Dynamics of a many-particle Landau-Zener model: inverse sweep 
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We consider dynamics of a slowly time-dependent Dicke model, which represents a many-body 
generalization of the Landau-Zener model. In particular, the model describes narrow Feshbach reso- 
nance passage in an ultracold gas of Fermi atoms. Adiabaticity is destroyed when a parameter crosses 
a critical value, even at very slow sweeping rates of a parameter. The dynamics crucially depends 
on direction of the sweep. We apply our recent analysis [A.P. Itin, P. Torma, arXiv:0901.4778vl] 
to the "inverse" sweep through the resonance, corresponding (in a context of Feshbach resonance 
passage) to dissociation of molecules. On a level of the mean-field approximation, the dynamics is 
equivalent to a molecular condensate formation from Bose atoms within a two-mode model. Map- 
ping the system to a Painleve equation allows us to calculate deviation from adiabaticity at very 
slow sweeps analytically. 



Destruction of adiabaticity is conceptually important 
issue in both quantum and classical mechanics. Nonadia- 
batic dynamics of many quantum single-particle systems 
can be effectively described in terms of the exactly solv- 
able Landau-Zener model (LZM), where a probabil- 
ity of transition from an initially occupied instantaneous 
ground state level to an excited one is exponentially small 
in a sweeping rate parameter. Numerous studies have 
been devoted to such type of models @. 

Many-particle quantum systems are intriguing objects 
that can exhibit nonlinear behavior and in some sense 
link quantum and classical worlds together. Interest to 
such systems has grown enormously due to progress in 
ultracold quantum gases experiments @ . Destruction of 
adiabaticity in semiclassical models of many-particle sys- 
tems under slow driving has been intensively investigated 
recently @, S !, 0, & @ using various methods. It has 
been found that exponential LZM-type behavior for tran- 
sition probabilities is often replaced with power-laws of 
sweeping rates. Theoretical methods used for derivation 
of these power-laws often involve a number of approxi- 
mations those accuracy is difficult to control. 

We considered recently [l(| a time-dependent Dicke 
model (see Eq. [2] below) using concepts of classical adi- 
abatic invariants. Technically, a method which natu- 
rally links classical mechanics and quantum dynamics of 
many-particle systems is a so-called Truncated Wigner 
Approximation (TWA) @, HE d- Briefly, within this 
approximation one considers an ensemble of classical tra- 
jectories whose initial conditions are distributed in phase 
space according to a Wigner transform of an initial quan- 
tum state. One therefore is left with the problem of 
nonadiabaticity of individual classical solutions. In clas- 
sical mechanics, it is known that adiabaticity is destroyed 
when a phase point crosses or comes close to a separatrix 
of an unperturbed system [l5j . There is a method for 
calculating changes in approximate adiabatic invariants: 
separatrix crossing theory [3, 0, EH • The method has 
a drawback: usually it is not applicable in case a separa- 
trix crossing happens close to or at the time of bifurca- 
tion, i.e. in case of very small initial action of a classical 
trajectory. In [10], we managed to overcome this draw- 



back by mapping the system close to a bifurcation to a 
Painleve equation, those asymptotics are known. A for- 
mula for change in the action we derived there allowed 
to describe mean values of number of bosons created at 
the resonance passage, and its distributions. In this Brief 
Report we consider the " inverse" sweep through the reso- 
nance in the same model. We also discuss other physical 
realizations of the semiclassical model, in particular as- 
sociation of Bose atoms in diatomic molecules within a 
simple two- model model @. 

Let us consider the time-dependent Dicke model writ- 
ten in the following form 



c-aepe 



H = -A(t)tfb + A(t)S z 



(1) 



where is the coupling strength, S = S x ± iS v are 

spin operators, w and b are creation and destruction op- 
erators of a bosonic mode, A(t) = ±ei is detuning, e is the 
sweeping rate of the bosonic mode energy. A(t) = +et 
for a "forward" sweep through resonance, and — et for a 
"backward", or "inverse", sweep. The spin value S is 
macroscopically large S = N/2 3> 1. With N = 2, one 
recovers the standard Landau-Zener model. We set g = 1 
for convenience in this paper. 

For a "forward" sweep through the resonance consid- 
ered in |10| . we started in the distant past with some 
small initial number of bosons Nt(t)\_ oc = (6^&)(t)|_ 00 = 
N- = ri-N, and calculated the final number of bosons 
Nb(t)\+oc = n f N by averaging over an ensemble of clas- 
sical trajectories. For a backward sweep, only numerical 
results were presented in Altland et al Q , consistent with 
the earlier work of Our method enables to explain 
them analytically with high accuracy. The classical limit 
of (TT]), obtained using the c— number formalism [7[,(9(, 
is 



H 



-771 - 



2nVl — ncoscj), 7 = ±2et, 



(2) 



where n corresponds to the number of bosons, and <f> ls 
the canonically conjugated phase. For a forward sweep, 
7 = 2ei, while for a backward sweep 7 = — 2et. Note 
that by means of a trivial change of variables [ijj], the 
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FIG. 1: Phase portraits of the Hamiltonian ([2]). From (a) to 
(f): 7 = -20, -4, -1.8, -0.9, 0., 0.9, 1.8, 20, 100 correspond- 
ingly. Saddle points are denoted by asterisks, the bold trajec- 
tory is the separatrix. Shaded areas illustrate definitions of 
the classical actions (see main text). 



Hamiltonian becomes the same as analyzed in [sj within 
the same context, and can be derived from semiclassical 
Bloch equations of [6]. A mean-field model of associa- 
tion of bosonic condensates into a molecular condensate 
considered in || is also equivalent to ((2|) , with n = 1 cor- 
responding to "all-atom" mode (i.e., on a mean-field level 
association of bosonic atoms is analogous to dissociation 
of molecular condensate of Fermi atoms ) . 

Let us recall the classical phase space of the Hamilto- 
nian ([2]) at fixed values of the parameter 7, as explained 
in [§]. If 7 < — 2, there is only one stable elliptic point 
on the phase portrait (at <j) = tt and n not far from 1, see 
Fig. la, b). At 7 = —2, a bifurcation takes place, and at 
I7I < 2 phase portraits look as shown in Figs.lc,d,e,f,g. 
There are two saddle points at n = 0, cos <f> = —7/2, and a 
newborn elliptic point at 4> — 0. The trajectory connect- 
ing these two saddles separates rotations and oscillating 
motions and it is called the separatrix of the frozen sys- 
tem [9(. Finally, at large positive 7, again there is only 
one elliptic stationary point at <f> — 0, and n close to 1. 
The classical action is defined as in 0, Q and is shown 
graphically in Figs la,e,i: shaded areas, divided by 27T. 
At 7 = —00, the action coincides with n: I — n. At 
7 = +00, we have I = 1 — n. For a forward sweep, 
starting at t = —00 with certain small n — n_, in the 
adiabatic limit we get n + = 1 — rt_ at t = +00. Analo- 
gously, for a backward sweep, starting from a n_ « 1, in 
the adiabatic limit we get n + = 1 — rt_ at t — +00. De- 
struction of adiabaticity in both cases is mainly caused 
by motion near the saddle points located at n — 0, i.e. 
no matter in what direction a sweep happens, region of 
the loss of adiabaticity is the same. For a forward sweep, 
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FIG. 2: Change in the classical action for the backward sweep 
with zero initial action (i.e., AI — If — J_ = If). Inset shows 
change in the action for small values of e and a linear fit 
which gives If/e = 0.220631. Theoretical prediction is 
^ w 0.220636. 



a classical phase point with J_ = (i.e., located initially 
at n = 0) always remains there. Therefore, when ap- 
proximating an initial quantum state by an ensemble of 
classical trajectories, one needs to carefully take into ac- 
count quantum fluctuations. For a ground state as an ini- 
tial state, initial actions of the classical ensemble ~ 1/N. 
For a backward sweep, the situation is different. A phase 
point with J_ = is not stationary and therefore in the 
semiclassical limit, one may merely study its classical dy- 
namics. For small sweeping rates, dynamics remains adi- 
abatic until the phase point approaches a saddle point 
at n = 0. Then, following [lfj, close to a saddle point 
the system can be mapped to a Painleve equation. To 
this end, one introduces new variables P — 2^/ncos<j)/2, 
Y = 2y / nsin(/>/2. After the transformations and approx- 
imations described in [l(| one obtains the Hamiltonian 
H = — — h ^ — \~^Ti where s is a new time variable 
and P,Y are canonically conjugated variables. One can 
see that the effective Hamiltonian leads to the Painleve-2 
equation [H, HI 



d 2 Y 



sY - 2Y 3 



(3) 



Asymptotics of the Painleve-2 equation were investigated 
by Its and Kapaev 1201 . At s — > — 00 the asymptotic 
solution to © is [2(yp] 



Y(s) = a(-s) * sin I ^(-s) 3/2 + -a 2 ln(-s) + ip 



and in the limit s — > +00 it is 

Y(s) = ±y|±p(2.)-i cos (^.^ - ^ 2 ln( S ) 

where (a,<p) and (p, 6) are the "action-angle" variables 
characterizing the solutions in the limits s — ► ±00. These 
variables are related by the following relations [2(J Hl[ : 



P 



I m . 1 + |P| 2 



7T 2|Im(p)| 
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6 = -- + -p 2 m2-argr( ip 2 )-arg(l+p 2 ), (4) 

P = (exp(7TQ! 2 )) X 

( 3 2l „ TT ^ / a 2 \ \ 

x exp I i— a ln2 — — i — I argl I i— J — up J , 

where arg(-) denotes the argument of a complex num- 
ber, and r(x) is the gamma function. As s — > ±oo, the 
adiabatic invariant I p of the system ([3]) approaches the 
quantities I~ or 7+ which are defined (to the lowest order 
terms) as 

Ip = l7-> Ip = (5) 

Relations © , © allow to calculate the change in the ac- 
tion both for a forward [1 01 ] and backward sweep. 

It turns out that the backward sweep through the res- 
onance has very different behaviour from the forward 
one. Most importantly, in the limit J_ — > change in 
the action becomes phase-independent. That is, at large 
J_ 3> e predictions of the Painleve mapping correspond 
to the separatrix crossing theory calculations. However, 
as 7_ decreases, the change in the action approaches the 
phase independent value, which is linearly proportional 
to e. 

Indeed, for a classical trajectory with I_ = from 
(SD,© it follows [HI 

/ p - = iln(l + b| 2 ) (6) 
For the backward sweep, I~ corresponds to //, while 



Ip corresponds to For the initial action /_ = we 
have \p\ = 0. Taking into account all the transformations 
leading to the Painleve equation, we obtain: 

AI = I f = ^e (7) 

TT 

This is the main result of the present paper. Numerical 
calculations (see figure 2) reproduce the coefficient — 
with 5-digit accuracy. In Ref . [5] , a linear power-law for 
deviation from adiabaticity was also found numerically, 
however the accuracy of the analytical solution reported 
was lower than in our case. 

To conclude, we found that the mapping of the many- 
particle Landau-Zencr model to the Painleve equation 
[Hj] enables to describe not only " forward" sweep through 
the resonance, but also the reverse driving process. The 
latter process arises, for example, when the model is used 
to describe formation of diatomic molecules from Bose 
atoms, or dissociation of diatomic molecules on Fermi 
atoms. Our results are relevant for accurately describing 
the Feshbach resonance passage in ultracold Fermi and 
Bose gases. Interferometric experiments with ultracold 
quantum gases [24j is another possible area of applica- 
tions. Finally, we believe the method will have important 
applications in the field of dynamics of quantum phase 
transitions (QPT) [H, ^ . 
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